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ABSTRACT 

Two consecutive transits of planetary companion OGLE-TR-lllb were observed in the I band. 
Combining these observations with data from the literature, we find that the timing of the transits 
cannot be explained by a constant period, and that the observed variations cannot be originated by 
the presence of a satellite. However, a perturbing planet with the mass of the Earth in an exterior 
orbit could explain the observations if the orbit of OGLE-TR-lllb is eccentric. We also show that 
the eccentricity needed to explain the observations is not ruled out by the radial velocity data found 
in the literature. 

Subject headings: planetary systems — stars: individual (OGLE-TR-111) 
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1. INTRODUCTION 

The observations of transiting extrasolar planets have 
produced some of the most interesting results in the 
study of other planetary systems. Their orbital con- 
figuration have permitted the first direct m easurements 
of radius, temperature, and composition (jSwain et alJ 
l2008t [Harrington et alj|2007l and references therein), all 
of which are critical to constraining the interior and ev o- 
lution models of extrasolar planets (e.g. iFortnev! 12008). 

It has been further realized that the presence of vari- 
ations in the timing of transits can be attributed to 
otherwis e undetectable planets in the system (see, for 
example. lMiralda-Escuddl2002t iHqlman fc Murravll2005t 



Ago! et al l 2005; H evl fe Gladma n 2007: Ford & H olmanl 



20071: ISimon et all I2007D . iDeeg et al.l (j2008D and 



Ribas et al.l ( 2008T ) reported indirect detections of un- 



seen companions by monitoring eclipse timing of the bi- 
nary stellar system CM Draconis (1.5 Mj to 0.1 M Q 
candidate) and variations in the orbital parameters of 
the planetary system around GJ 436 (5 M$ compan- 
ion), respectively. Howe ver, this last case ha s been re- 
cently argued against bv lAlonso et alJ (I2008D. Besides , 
recently-discovered transiting planets (jPont et alJl2007t 
lUdalski et alJ |2008[ ) exhibiting shifts in their radial ve- 
locities are promising new candidates to search for varia- 
tions in the tim i ng of their transits. On the other hand, 
ISteffen fc Agoll (|2005l ) found no evidence of variations 
in the timing of transits of the TrES-1 system, after 
analysing data for 12 transits. Also, after monito r ing 15 
transits of the star HD 209458. iMiller-Ricci etafl (|2008l ) 
were able to set tight limits to a second planet in the 
system. 

Here we report a significant detection of variability in 
the timin g of the trans its of extrasolar plane t OGLE- 
TR-lllb (lUdalski et aHl200l iPont et al][200l and dis- 
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cuss its possible causes, including a second unseen planet 
OGLE-TR-lllc. 

In a previous work (jMinniti et al. 2007) we reported 
a single transit observed in the V band which occurred 
around 5 minutes be fore the expected t ime obtained us- 
ing the ephemeris of I Winn et al.l (|2007l . hereafter W07), 
but the result was inconclusive since it had a 2.6-cr sig- 
nificance. In the present work we analyse data of two 
consecutive follow-up transits of the same planet. 

Section[2]presents the new data and the reduction pro- 
cedures, in Section [3] we describe the technique used to 
measure the central times of the transits. Finally, in 
Section [4] we present our results and discuss their impli- 
cations. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed two consecutive transits of planetary com- 
panion OGLE-TR-lllb in the I band with the FORS1 in- 
strument at the European Southern Observatory (ESO) 
Very Large Telescope (VLT). The observations were ac- 
quired during a Director's Discretionary Time run during 
the nights of December 19 and December 23, 2006. Since 
the orbital period of OGLE-TR-lllb (P = 4.01444 days) 
is almost an exact multiple of Earth's rotational period, 
those were the last events visible from the ESO facilities 
in Chile until May 2008. 

FORS1 is a visual focal-reducer imager who had a 
2048x2048 Tektronik CCD detector and a pixel scale of 
0.2 arcsec/pix. For the observations, a nearby bright star 
was moved outside the field of view, leaving OGLE-TR- 
111 near the center of the north-eastern quadrant. The 
chosen integration time of 6 seconds was the maximum 
possible to avoid saturation of the star in case of excel- 
lent seeing. A total of over 9 hours of observations were 
obtained during the second half of both nights. Dur- 
ing the first night the seeing remained stable below 0.6" , 
but it oscillated between 0.6" and 1.4" during the second 
night. Observations finished near local sunrise produc- 
ing a non-centered bracketing of the events and an addi- 
tional source of scatter as the sky background increased 
near sunrise. 

We used the ISIS package (jAlard fc Luptonl 119981 : 
lAlardl[2000l ) to compute precise differential photometry 
with respect to a reference image in a 400 x 400 pix sub- 
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TABLE 1 

Orbital and physical parameters for system OGLE-TR-111. 



Time from miill nmsil lUiysj 

Fig. 1. — Relative flux during two consecutive transits of plan- 
etary companion OGLE-TR-lllb. Except for those mentioned in 
the text, no points were discarded. In the upper (lower) panel we 
present data taken on the night of December 19 (23) 2006. The 
residuals with the error bars are also shown. The dashed line rep- 
resents the displaced zero for the residuals, and the (red) solid line 
is the best fit model. Note how the errors increase at the end of the 
second night due to the increase of the background noise caused 
by dawn. 

frame. The reference image was obtained combining the 
10 images with best seeing, which produced an image 
with FWHM « 0.46 arcsec. The resulting subtracted 
images were checked for abnormally large deviations or 
means significantly different from zero; an image from 
the first night and three images from the end of the sec- 
ond night were discarded in this way, leaving a total of 
488 images. 

Aperture photometry was performed o n the differenc e 
images using IRAF DAOPHOT package (|Stetsodll987h . 
which was found to give better results than the ISI S pho- 
tometry routine phot. csh (se e lHartman et al.ll2004h . In 
agreement with iGillon et al.l (|2007f ) . we found that the 
scatter increased rapidly with aperture size, although in 
our case the transit amplitude remained constant (within 
a 0.1% level). We therefore choose a 5-pixels aperture 
since our goal is to obtain precise measurements of the 
central times of transits, and therefore the relevance of 
obtaining the correct amplitude is diminished. 

The uncertainty in the difference flux was es- 
timated from the magnitude error obtained from 
DAOPHOT/APPHOT, which uses Poisson statistics, 
and considers the deviation in the sky background. The 
flux in the reference image was measured using PSF- 
fitting photometry with DAOPHOT/ALLSTARS. The 
systematic error introduced by this measurement is stud- 
ied further in Sect. [3] 

To remove possible systematics effects from the light 
curves we employed the signal reconstruction meth od of 
the Trend Filtering Algorithm (jKovacs et al.ll2005h . We 
refer readers to this paper for a description of the method 
as well as for an illuminating discussion of the possible 
causes of systematics effects. We chose light curves of 19 
stars distributed as uniformly as possible around OGLE- 
TR-111 as template light curves, and checked them for 
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Fig. 2. — Probability density distributions for the central times 
of the transits obtained from the MCMC simulations. The thick 
vertical solid line indicates the median of the distribution, and the 
dotted lines mark the upper and lower 68% confidence limits. The 
solid (blue) curve is a Gaussian distribution with the same mean 
and standard deviation as the data. 

obvious variability or uncommonly large scatter. The al- 
gorithm was iterated until the relative difference in the 
curves obtained in two successive steps was less than 
10~ 5 . The resulting science light curves for both nights 
are shown in Fig. Q] The standard deviation before the 
transit of the second night is 2.65 mmag, almost reaching 
the photon noise limit of 2.55 mmag. 

3. MEASUREMENTS 

Planetary and orbital parameters, including the cen- 
tral times of transits, were fitted to the OGLE-TR-111 
light curve. The model used consisted on a perfectly 
opaque spherical planet of radius R p and mass M p , or- 
biting a limb-darkened star of radius R s and mass M s 
(|Mandel fc Agoll |2002) in a circular orbit of period P 
and inclination i. We considered a quadratic mo del for 
the lim b-darkening, with coefficients taken from I Claret! 
(2000) for a star with T cff = 5000 K, \ogg = 4.5 cms" 2 
and [Fe/H] = 0.2 and microturbulent velocity £ = 2 
km/s. The mass of the pl anet and the s tar w ere fixed 
to the values reported by ISantos et al.1 (|2006h . M s = 
0.81 Mq and M p = 0.52 Mj up . The remaining five pa- 
rameters for the model: R p , R s , i and the central time 
of each transit (T c i and T C 2) were adjusted using the 488 
data points of the light curve. 

The parameters were obtained by minimizing the 
X 2 statistic us i ng th e downhill simplex algorithm 
(|Nelder fc Meadl 119651 ) implemented in the Scipy li- 
brary 6 . We present the parameters in Table [lj and the 
best-fit model and the residuals in Fig. [TJ Note that, ex- 
cept for the planetary radius and the time between first 
and last contact, the parameters reported in Tabled] are 
in agreement with previously published values (see Sect. 

a. 

The uncertainties in the parameters were estimated us- 
ing the Markov Cha in Monte Carlo metho d, whic h is de - 
scribed in detail by iTegmark et all (|2004l ). iFordl (|2005l ) 

6 http://www.scipy.org 
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and iHolman et alj (|2006l ) . We constructed chains with 
500.000 points each, and discarded the first 100.000 to 
guarantee convergence. The jump function employed was 
the addition of a Gaussian random number to each pa- 
rameter, and a global scaling of the sigma of the random 
Gaussian perturbations was adjusted after convergence 
was reached so that between 20% and 30% of the jumps 
were executed. 

In this manner, we built five independent chains and 
found that the mean values and the confidence inter- 
vals of the parameters (computed as described below) 
are in excellent agreement for all chains, a sign of good 
convergence. Besides, the correlation length, defined as 
the number of steps over which the correlation function 
fsee lTegmark et aLll2004 Appendix A) drops to 0.5 was 
about 80 for the central times of the transits, and around 
800 for the high ly cov ariant parameters R p , R s and i, in 
agreement with IW0"7l . This produces an effective length 
of about 5000 for T c \ and T c2 , a sign of good mixing. 

For each chain we took a random subset of 5000 val- 
ues (the effective length) of the central times and test 
the hypothesis that the sets were drawn from identi- 
cal populations using the Wilcoxon's rank sum test (see 
iFrodesen etHl 119791 §14.6.9). For all cases the test 
statistic (which is approximately Gaussian) falls within 
2.5-sigma of the expected value, and therefore the hy- 
pothesis cannot be discarded for significance levels below 
« 1.2%. 

Fig. [2] shows two representative probability density dis- 
tributions corresponding to the two central transit times 
and Table [T] reports the median and the upper and lower 
68% confidence limits, defined in such a way that the 
cumulative probability below (above) the lower (upper) 
confidence limit is 16%. As a solid curve we plot the 
Gaussian probability density having the same mean and 
standard deviation as the data. 

To test the robustness of our results, the fit was re- 
pe ated fi xing the values of R p , R s and i to those reported 
bv IW071 (Rp = 1.067 R Jup , R s = 0.831 R Q , i = 88.1 
degrees) and including the out-of-transit flux as an ad- 
justable parameter. The obtained times for the center of 
the transits are in agreement with those reported above. 
The sa me res ults are obtained if only R s is fixed to the 
value of lW07l 

Additionally, to check that the systematics-removal 
procedure does not modify the shape of the light curves, 
we also measured the central times in the original curves 
obtained with aperture photometry. Again, the obtained 
values are in excelent agreement with the ones presented 
above, and the errors computed with MCMC are larger 
by a factor between 1.04 and 1.99, depending on the pa- 
rameter, as expected. 

Possible systematic errors may be introduced by the 
choice of the stellar mass, the orbital period — which af- 
fects the determination of the orbital radius — , the model 
for the limb darkening, and the flux in the reference im- 
age. To study these effectes we obtained new fits to 
the data varying the fixed parameters and the function 
for the limb darkening. The stellar mass was varied by 
±10%, the photometry in the reference image was varied 
by ±0.1 mag and the orbital period by ±10 a (see Eq.|5J). 
The coefficients for the quadratic limb-darkenning model 
w ere adjusted fr om the data instead of fixed to the values 
of IClaretl (|2000T ) and, additionally, a linear limb darken- 
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Fig. 3. — Observed-minus-calculated times (in minutes) for the 
transits of planet OGLE-TR-lllb in front of its host star. The 
filled circles are the n ew transits presented in this work, the empty 
circles arc from W07 and the empty square is the transit presented 
bv lMinniti et alT (:2007), which has been reprocessed for this work. 

ning model was considered, bo th fixing the linear coeffi- 
cient to the value computed bv lClaretl (|2000f ) and adjust- 
ing it as part of the fit. In all cases, the variation in the 
central times of transit was smaller than the uncertain- 
ties reported in Table [TJ We therefore conclude that the 
values obtained for the central transit times are robust. 

4. RESULTS AND DISCUSSION 

We fitted a straight line to the centra l ti mes of the two 
transi ts together with those from lW07l and lMinniti et all 
(2007). The central time of this last transit (T c< vimos) 
has been remeasured using the procedure described 
above and the result is shown in Table [T] In this way 
we obtained a new ephemeris for the transit times: 

T c = 2454092.80607 ± 0.00029 [HJD] (1) 
P = 4.0144540 ±0.0000038 days , (2) 

with correlation coefficient p = 0.785. The reduced \ 2 is 
9.04, indicating a poor fit. Note that the v alue of the pe- 
riod is consistent with the value reported by |W07l . The fit 
was repeated including a point for the OGLE data, and 
we also obtained the period from a simu l taneous fit to 
all the available photometry (OGLE. rW07l . lMinniti et al.l 
(|2007f ) and this work). In both cases the obtained value 
is in excellent agreement with the one reported above. 

In Fig. [3] we plot the residuals of the fit. It is clear 
that the observed-minus-computed (O-C) values are not 
consistent with a constant p eriod s ince the VIMOS tran- 
sit, one of the transits from IW0"7l and one of the FORS 
transits lie -3.29-cr, 2.79-cr and -2.52-ct away from zero, 
respectively. However, the data available to date are not 
enough to determine the nature of these variations. Nev- 
ertheless, we have been able to discard a few possibilities 
and study some others. We present some preliminary re- 
sults here and defer a more detailed study for a future 
work. 

First, the hypothesis of an exomoon seems unlikely, 
since the mass needed to produce the observed O-C am- 
plitude is at least one twenty-sixth of the planetary mass 
if the moon is at a Hill radius from the planet. How- 
ever, at this distance the moon system is expected to 
be unstable. For moons closer to the planet, the needed 
mass increases. These are extreme values when compared 
with the S olar System , where this ratio never exceeds 
2.5 x 10- 4 (lCoxH2000lh 

On the other hand, several planetary system config- 
urations reproduce the observed trend. The equations 
of motion for the three-body problem were solved with 



4 




Date [JD-2453000] 



Fig. 4. — Radial velocity measurements from Pont ct al. II2004D 
together with the best fit (solid line), and the corresponding ±1<t 
curves (dotted lines). Also shown is the fit for e = (dashed line). 

the Bulirsch-S toer algorithm implemented in the Mer- 
cury package (|Chamberslll999f ) using different sets of or- 
bital parameters for the perturbing planet, and the re- 
sults were compared with the observations. A particu- 
larly interesting solution is that an exterior Earth-mass 
planet near the 4:1 resonance produces the observed am- 
plitude and periodicity in the O-C times, if the orbit of 
TRlllb is eccentric (e = 0.3). On the other hand, the 
mass of the perturber planet must be at least around 4 
Mj up if the orbit of the interior planet is nearly circular. 
This shows the importance of accurately measuring the 
ecentricity of the interior planet through R V data or mea- 
surem ents of the planet occultation (see iDeming et al.l 
l2007h . 



In the discovery paper by iPont et al.l (|2004D . the or- 
bital solution was obtained by fixing the eccentricity of 
TRlllb to zero. Although this is reasonable for a single 
planet in a close orbit to the star, since circulariza tion is 
very effective in those conditions (see, for example, |ZahnJ 
Il977f ). a second planet can perturb the orbit of the first 
one, increasing its eccentricity. Th erefore, we reanaly sed 
the radial velocity (RV) data from lPont et alj (|2004h . in 
order to constrain the possible eccentricity of the system. 
We found that the data are compatible with an eccentric- 
ity of 0.3, with a reduced \ 2 of about 0.4 (for 5 degrees 
of freedom, see Fig. rjj compared to the value of 0.7 for 
a circular orbit, as reported in the original paper. 

Additionally, note that the 1.55-cr difference between 
the transit Jength presented in Table [1] and that re- 
ported by IW07I might indicate a ch ange in the incli- 
nation angle of OGL E-TR-lllb (see iRibas etall 120081 : 
lMiralda-Escudel[2002f ) which could in principle help con- 
strain the parameters of the perturber planet. 

Future observations are warranted in order to pinpoint 
the origin of the variation in the period of this interesting 
planet. 
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